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Abstract. We show that, under suitable assumptions, Poincare recurrences of a 
dynamical system determine its topology in phase space. Therefore, dynamical systems 
with the same recurrences are topologically equivalent. 

Introduction: From the beginning of what is today called dynamical systems theory, 
there has been great interest in the following problem: what can be inferred about a 
system if only partial information about its state at each moment in time is given? This 
matter is highly relevant in physics, where often not all state variables of a system but 
only one observable can be measured. To make the following arguments more precise 
we introduce the following definitions. 

Definition 1.1. Let M be a compact manifold. A dynamical system on M is a diffeo- 
morphism (ft : M — ► M (discrete time) or a vector field X on M (continuous time). In 
both cases the time evolution corresponding to an initial position x e M is denoted by 
(j)t(xo). In the discrete case t G N and (pi = ((f)) 1 . In the case of continuous i £ R and 
t — ► (f>t(xo) is the X integral curve through Xq. 

Definition 1.2 : An observable is a smooth function y : M — > R. 

In the late 70's it was known that a procedure of delay embedding could be used to 
reconstruct the "Geometry [of an attractor] from a Time Series" [I]. This reconstruction 
could then be used to estimate dynamical invariants such as Lyapunov exponents. In 
1981, Takens proved three theorems [2], the first of which reads (using definitions 1.1. 
and 1.2.) 

Theorem 1.1: Let M be a compact manifold of dimension m. For pairs 

((j),y),4> : M — > M a smooth diffeomorphism and y : M — ► M a smooth function, it 

is a generic property that the map $($, y ) '■ M — > R 2m+1 , defined by 

^ y (x) = (y(x),y( ( p(x)),..., y (<p 2m (x))) (1) 

is an embedding; here "smooth" means at least C 2 . 

As mentioned above, the delay embedding procedure (pQ) lies at the heart of time se- 
ries analysis and several authors have extended Takens' Theorem [HI 0]. These theorems 
are virtually always used when nonlinear time series analysis is applied to experimental 
data p]. As the reconstructed attractor is diffeomorphic to the original attractor, it 
can be considered to be the same system after a smooth co-ordinate change. Such a 
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change leaves dynamical invariants, such as entropies, fractal dimensions and Lyapunov 
exponents unchanged. 

In this paper we study a similar situation, in which it is not an observable in the 
sense of definition 1.2. which is known, but rather the Poincare recurrences of the 
system. The theorems we report are not primarily intended for applications in time 
series analysis, but rather to stress the importance of Poincare recurrences. However, 
there are some applications, e.g. in the field of data analysis namely for the method 
of recurrence plots |6j or for the reconstruction of protein structures [7] for which the 
results are relevant. 



Recurrence is a fundamental characteristic of many dynamical systems. The 
concept of recurrence was introduced by Poincare in a seminal work in 1890 [9J, for 
which he won a prize awarded by King Oscar II of Sweden and Norway. Therein, 
Poincare not only described the "homoclinic tangle" which lies at the root of the chaotic 
behaviour of orbits, but he also introduced (as a by-product) the concept of recurrences 
in conservative systems. When speaking about the restricted three body problem he 
reported :" In this case, neglecting some exceptional trajectories, the occurrence of which 
is infinitely improbable, one can show, that the system recurs infinitely many times as 
close as one wishes to its initial state." In a larger context recurrences make part of one 
of three broad classes of asymptotic invariants |8J: (i) growth of the number of orbits of 
various kinds and of the complexity of orbit families^, (ii) different types of recurrences, 
and (iii) asymptotic distribution and statistical behaviour of orbits. The first two classes 
are of purely topological nature. The last class is naturally related to ergodic theory. 
Of the different types of recurrences which form part of the second class of invariants, 
the Poincare recurrence is of particular interest to this work. It is based on the Poincare 
Recurrence Theorem [8] (Theorem 4.1.19): 

Theorem 2.1. Let T be a measure-preserving transformation of a probability space 
(X, jj) and let A C X be a measurable se|§|. Then for any natural number N G N 



Due to the great relevance of this theorem for this paper we give here a rather short 
(standard) proof of this theorem: 

Proof: Replacing T by T N in Eq. ([2]) , we find that it is enough to prove the statement 
| An important invariant of the orbit growth is the topological entropy 

§ Here ^ is a Borel measure on a separable metrisable space X. Note that these assumptions are 
rather weak from a practical point of view. Such a measure preserving function is obviously given in 
Hamiltonian systems and also for all points on (the w-limit set of) a chaotic attractor. 
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for N = 1. The set 

A:={xeA \ {T n (x)} nm c X\A} = An ^f] T- n (X\A)\ 
is measurable. Note that 

T- n (A)nA = ® 

for every n, because if we suppose, on the contrary, that T~ n (A) H A = B with 5^0, 
this implies T n (B) C A. This is inconsistent with the definition of A, because B C A. 
Also note that 

r n (i)nr m (i)=i Vm,fieN, 

because if we assume, on the contrary, that T~ n (A) H T~ m (A) = B with 5^0, this 
implies T n (B) C *4 and T m (B) C ^4. Without loss of generality, we assume that m > n. 
Then r(B) = C C i and T m (S) = T m ~ n (T n (B)) = T m ~ n (C) C „4, which is again 
inconsistent with the definition of A. 

Furthermore, fi(T~ n (A)) = ^(A) since T preserves \i. Thus = because 

too \ oo oo 

|J T-"\A) = ^ M (t-(^)) = £ □ 
n=0 / n=0 n=0 

That means, that if we have a measure preserving transformation, the trajectory will 
eventually come back to the neighbourhood of any former point with probability one. 
However, the theorem only guarantees the existence of recurrence but gives no indication 
of how long it takes the system to recur. Especially for high dimensional complex systems 
the recurrence time might be extremely long. 

The first return of a set is defined as follows: if A C X is a measurable set of a measurable 
(probability) dynamical system {X,fi,T}, the first return of the set A is given by 

t(A) = min {n > : T n A n A ^ 0} . 

Generically, for hyperbolic systems the recurrence or first return time appears to exhibit 
certain universal properties [TO] : 

(i) the recurrence time has an exponential limit distribution; 

(ii) successive recurrence times are independently distributed; 

(iii) as a consequence of (0} and (EJ) , the sequence of successive recurrence times has a 
limit law that is a Poisson distribution. 

These properties, which are also well-known characteristics of certain stochastic systems, 
such as finite aperiodic Markov chains [TTJ, [T2l H3] , have been rigorously established for 
deterministic dynamical systems exhibiting sufficiently strong mixing [HI CDS [16]. ^ ne y 
have also been shown valid for a wider class of systems which have to be hyperbolic [17]. 
Recently, recurrences and return times have been studied with respect to their statistics 
pEEICEH] and have been linked to various other basic characteristics of dynamical systems, 
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such as the Pesin's dimension [20], the point- wise and local dimensions [21], [22], [23] or 
the Hausdorff dimension [24]. Also the multi-fractal properties of return time statistics 
have been studied [25] [26]. Furthermore, it has been shown that recurrences are related 
to Lyapunov exponents and to various entropies [27] [28]. They have been linked to rates 
of mixing [29], and the relationship between the return time statistics of continuous and 
discrete systems has been investigated [30] . 

An alternative and very practical approach has been introduced by Eckmann et al. in 
1987 [31]. Instead of studying the (first) return time statistics, they introduced the 
recurrence matrix 

# M - = e( £ Hl^o)-^o)ll), (3) 

where 4>i(xo), (fij(xo) are the points on the continuous or discrete orbit of a dynamical 
system at discrete times i and j, i.e. Eckmann et al. considered vector valued time series. 
If the (Euclidean) distance between these two points is smaller than the (arbitrary but 
small) threshold distance e, the Heaviside function G(-) leads to R it j = 1. If their 
distance is larger than e we have Rij = 0. Note that we will adapt the definition of 
the Heaviside function for which 6(0) = 0. In principle, the definition 6(0) = 1 would 
also be possible. Based on these definitions, R it j is a binary matrix, which contains 
information about the recurrences of the dynamical system. 

This means that the orbit of the dynamical system is mapped to a binary matrix, 
which indicates that, in a sense, much of the information these orbits contained is lost. 
However, the recurrence matrix and the derived recurrence plots, which are a graphical 
representation of the recurrence matrix, have been applied very successfully to analyse 
time series. From the matrix many dynamical invariants, such as generalised Renyi 
entropies and dimensions and also the mutual information, could be estimated [32j. 
Recently, it has been shown that for one-dimensional systems the time series can be 
reconstructed from the recurrence matrix up to a smooth observation function [33j. 
Later a numerical algorithm was used to reconstruct the attractor from the recurrence 
matrix of various systems, e.g. the Lorenz attractor [34]. An interesting feature of these 
algorithms is that they are based on purely geometrical considerations, which is different 
from delay embedding procedure in Takens' Theorem. 

These results suggested that once the recurrences of a system are known (in the form 
of the recurrence matrix in the limit of long time series) its dynamics are determined. 
In the next section we will prove a theorem, which will show that under rather weak 
assumptions it is possible to "reconstruct the geometry" of an attractor if the distance 
inequalities are known for each pair of points. 

1. Proof of the Reconstruction Theorem 

To investigate how much recurrences tell us about the dynamics of a system, we suppose 
that we have a subset M of M n for some n, and we know for each pair of points their 
distance inequality, i.e. we know if d(x, y) < e, or if d(x, y) > e. In terms of dynamical 
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systems theory, these might be the points on the attractor in phase space, together with 
their "recurrence matrix". Suppose we then map all points of M onto another set N 
by a function /, such that the resulting distance inequalities are the same, i.e. we have 
d(f(x),f(y)) < e if and only if d(x,y) < e. Now we have two sets (with an induced 
topology inherited from R ra ), which have the same distance inequalities, or from the 
point of view of dynamical systems theory, the same recurrences. In the context of the 
dynamical considerations, we impose a fairly natural additional condition: that any two 
elements of M can be distinguished by comparing their distances (relative to e) from 
enough elements of M. More precisely, if the "separation" condition of the Theorem 
(and its Corollary) below fails for M, then there are distinct points x,y G M such 
that the closed ball (in M) of radius e with centre x is contained in the corresponding 
closed ball in M with centre y. In the following proofs, we need to distinguish carefully 
between strict inequalities and cases where equality might be allowed. In the practical 
applications, this distinction seems unlikely to matter. 

A corollary, which is a consequence of the following theorem will show, that / must 
then induce a homeomorphism. 

RECONSTRUCTION THEOREM: Let M, N be closed subsets of W l for some 
n < oo. Let f : M — > iV be a surjection with the property that for some fixed e > 0, 
we have d(f(x),f(y)) < e if and only if d(x,y) < e, where d is the usual metric. 
Suppose also that, for each x 7^ y G M, there is an element z G M with d(x, z) > e 
and d(y, z) < e, and similarly, whenever u 7^ v G N, there is an element w & N with 
d(u, w) > e and d(v, w) < e. Then f is a homeomorphism. 

PROOF: Notice that / is injective, since for x,y,z above, d(f(x),f(z)) > e and 
d(f(y),f( z )) < £, so f(x) and f(y) are distinct. Hence, by symmetry, it suffices to 
prove that / is continuous. 

Notice, then, that for each x G M, 

M\{x}= [J B(z,e), 

{z£M:d(x,z)>e} 

where B(z, e) denotes the open ball in M with centre z and radius e, since x is outside 
this union, but, by hypothesis, any y different from x is in the union (it is at distance 
less than e from some z such that d(x, z) > e). 
Also, 

N\{f(x)}= |J B'(f(z),e), 

{zeM:d(x,z)>e} 

where B'(v, e) denotes the open ball in iV with centre v and radius e. To see this, we 
apply / to each side of the equality 

M\{x}= [J B(z,e) 

{z£M:d(x,z)>e} 
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and note that 

f(B(z,e)) = B'(f(z),e) 

by the hypotheses. 

Take a sequence (xk) of elements of M having limit x G M. We will prove that 
{f{ x k)) has limit f(x). Let K = B'(f(x),e) (since N is closed in R n , the closures in 
N and M. n coincide). Notice that K is closed and bounded, so is compact. For large 
enough k, we have Xk G B(x,s) and f(xk) G B'(f(x),e), using the hypotheses for the 
latter inclusion. 

Choose 5 > 0. Now 

N = B'(f(x),5)U |J #(/(*),£), 

{zeM:(i(a;,z)>£} 

and 

K = (K n 5)) u |J (ifn 

{2£M:ii(i,z)>£} 

a cover of K by open sets (in K). Since K is compact, there is a finite subcover, which 
must include (K fl B'(f(x), 5)), as none of the open sets B'(f(z), e) (where d(x, z) > e) 
contains f(x), (for d(x, z) > e yields d(f(x), f(z)) > e by hypothesis). 

Hence there is an integer s and elements v\, . . . , v s G M such that d(x, Vj) > e and 
d(f(x), f(vj) > e (for each j) and 

s 

K = (K f] B'(f(x), 5)) u\J(Kn B'(f( Vj ),e)). 

i=i 

Let 

a = mm{d(x, Vj) — e : 1 < j < s} > 0. 

For all large enough m, we have x m G B(x,e), f(x m ) G B'(f(x),e) and d(x,x m ) < a. 
Then (for each j) 

d(x, vj) < d(x, x m ) + d(x m , vj) < a + d(x m , Vj), 

so that 

d(x m , Vj) > d(x, Vj) — a > d(x, Vj) — (d(x, Vj) — e) = e, 
as a < d(x, Vj) — e. Hence d(f(x m ), f(vj) > e for each j. 
Since 

s 

K=(Kn B'(f(x),5)) u\J(Kn B'(f(vj),e)) 

j'=i 

and d(f(x m ), f{vj)) > e for each j, it must be the case that d(f(x), f(x m )) < 5. Since 5 
is arbitrarily chosen, (f(xk)) has limit f(x). Thus / is continuous, as required. □ 
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The following useful corollary is immediate, since compact subsets of M. n are closed, 
and because the separation property required in N follows automatically from the prop- 
erties of / and the fact that the corresponding property holds in M. 

Corollary 1: Let M, N be compact subsets ofM. n for some n < oo. Let f : M — ► iV 

be a surjection with the property that for some fixed e > 0, we have d(f(x), f(y)) < e if 
and only if d(x,y) < e and d(f(x),f(y)) > e if and only if d(x,y) > e, where d is the 
usual metric. Suppose also that, for each x ^ y G M , there is an element z e M with 
d(x, z) > e and d(y, z) < e. Then f is a homeomorphism. 

Interpretation 

The theorem and the corollary given in the last section show that distance 
inequalities restrict the possible structure of a set up to homeomorphic transformations. 
The recurrence matrix also contains information on the time evolution of the system. 
The indices i and j in Equation j3]) correspond to the times when the points </>j(xo) 
and 4>j(xo) are close to each other. This means that the recurrence matrix contains the 
information of whether or not the system at time j has recurred to the state it was in at i. 
If the system lives in its phase space M, the theorem and the corollaries guarantee that if 
we take all the recurrences into account, we can reconstruct a homeomorphic set N, once 
points with identical neighbourhoods are identified. We furthermore know the sequence 
in time in which the different points in the phase spaces are visited. We therefore 
can reconstruct a topologically conjugate dynamical system from the knowledge of the 
recurrences. This means that the reconstructed system is dynamically conjugate to the 
original one, in the sense of [35]. Both systems can then be considered to be identical 
up to a continuous change of the coordinate system. 

Conclusions 

Several conclusions can be drawn from these considerations. First of all, the 
theorem shows that under some assumptions recurrences determine a system up to 
topological conjugacy. Two systems with the same recurrence matrix are dynamically 
equivalent and can be considered to be linked by a continuous change of coordinates. 
This result also sheds light on standard techniques which are used to detect generalised 
synchronisation, where one basically considers simultaneous recurrences in both systems 
[36, 37]. These methods quantify the similarity of both recurrence matrices for small 
thresholds e, yielding an index for generalised synchronisation. Furthermore, the main 
result of this paper provides a theoretical basis for the numerical algorithm used in 
|34j to reconstruct the trajectory of the underlying system from its recurrence matrix. 
This reconstruction was shown to be possible for a rather large interval of values of the 
threshold e. In the case of the Bernoulli map, for example, the e-ball could cover up to 
half of the phase space. Also for continuous systems like the Rossler or the Lorenz system 
this reconstruction could be performed for ^-neighbourhoods covering a large fraction 
of the phase space. The quality of the reconstruction was practically independent of 
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e then, as even spatial structures much smaller than e could be reconstructed for long 
time series. 

The theorem proved in this paper also supports the method proposed in [38J to generate 
twin surrogate data from the recurrence matrix. 

Furthermore, our considerations suggest that it be might possible to show that / 
together with its inverse f~ l has to be differentiable, i.e. at least C 1 , and therefore 
induces a diffeomorphism between M and N . 

This work illustrates that, in a sense, recurrences provide an alternative description of 
a dynamical system. 
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